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Abstract.  Fundamental  understanding  of  dynamic  behavior  of  polycrystalline  ceramics  is  advanced 
through  constitutive  theory  development  and  computational  modeling.  At  the  mesoscale,  micro¬ 
structures  of  silicon  carbide  grains  (hexagonal  crystal  structure)  or  aluminum  oxynitride  grains  (cubic 
crystal  structure)  are  subjected  to  compression  or  shear  at  high  rates  with  varying  confining  pressure. 
Each  grain  is  resolved  by  numerous  three-dimensional  finite  elements,  and  behavior  of  each  grain  is 
modeled  using  nonlinear  anisotropic  elasticity.  Cohesive  fracture  models  and  post-fracture  contact  are 
included.  Normal  and  Weibull  failure  statistics  from  many  simulations  are  collected  and  analyzed. 
Results  demonstrate  effects  of  load  directionality,  confinement,  dilatation,  elastic  anisotropy  and 
elastic  nonlinearity,  and  grain  boundary  fracture  properties  on  macroscopic  (average)  failure  stresses 
for  loading  conditions  in  the  ballistic  regime.  Predictions  demonstrate  reasonable  agreement  with  data 
from  macroscopic  plate  impact,  unconfined  compression,  and  flexure  experiments. 

Keywords:  mesoscale,  ceramics,  silicon  carbide,  aluminum  oxynitride,  shock  compression,  fracture. 
PACS:  62.20.-x,  62.20.D-,  62.20.-mj,  81.05.Je. 


INTRODUCTION 

Dynamic  deformation  and  failure  behaviors  of 
high-strength  ceramics  are  of  interest  for  military 
and  industrial  applications.  Properties  of 
importance  in  the  ballistic  loading  regime  include 
hardness,  elastic  stiffness,  fracture  toughness, 
unconfined  compressive  strength,  dynamic  shear 
strength,  and  failure  probabilities  [1-6].  The 
present  work  investigates  two  ceramic  materials: 
silicon  carbide  (SiC)  and  aluminum  oxynitride 
(AlON).  Mesoscale  models  are  used,  wherein  each 
grain  in  a  polycrystalline  aggregate  is  resolved 
explicitly  in  a  numerical  simulation. 

This  paper  is  organized  as  follows.  First, 
models  for  single  crystals  and  interfacial  fracture 
are  listed.  Next,  three-dimensional  microstructures 
depicted  using  finite  elements  are  described. 
Results  of  simulations  considering  uniaxial  strain 


compression,  uniaxial  stress  compression,  and 
shear  boundary  conditions  are  reported.  Finally, 
statistics  of  dynamic  failure  are  analyzed. 

THEORY 

Essential  aspects  of  the  constitutive  theory  are 
given  here;  a  comprehensive  description  will  be 
published  later  [7],  Consider  hyperelastic  behavior 
of  single  crystals  [8-10]  within  a  given  polycrystal. 
Spatial  and  material  coordinates  are  related  by 

xa=xa(XA,t),  (1) 

where  t  denotes  time.  The  deformation  gradient  is 
FaA  =  dxa/dXA.  (2) 

The  symmetric  Lagrangian  elastic  strain  is 

EAB=i(FaAFaB-GAB),  (3) 

with  Gab  the  covariant  referential  metric  tensor. 


Shock  Compression  of  Condensed  Matter  -2011 
AIP  Conf.  Proc.  1426,  1039-1044  (2012);  doi:  10.1063/1.3686456 
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Assuming  isentropic  behavior,  internal  energy 
density  per  unit  reference  volume  U  is,  to  third 
order  in  strain, 

[j  =  lrABCDF  F  + 1  pABCDEF  p  p  p 

In  Voigt’s  notation,  second-  and  third-order  elastic 
constants  CABCD  <H>  Cafs  and  CABCDEF  <->  CaP/  . 

The  first  Piola-Kirchhoff  stress  is 

PaA  =JF~lA(7ab  =dU/dF,. 

—  Fa(C'ABCDF  I  j_  abcdef  t?  f  \  v  ; 

r .B  v.  ^ CD  '  2  ^  ^CD^EF/’ 

where  J  =  p0/  p  is  the  Jacobian  determinant  of  (2) 
and  <7ab  is  the  Cauchy  stress.  In  the  absence  of 
measurements  of  third-order  elastic  constants, 
pressure  derivatives  of  second-order  coefficients 
provide  the  following  approximation  [1 1,  12]: 

/~i  ABCDEF  B  /  d  ABCD  EF  .  d  /^iCDEF/^iAB 

C  ““T  (^C  G  G 

(6) 

+  jLCefabGcd\ 

where  p  =  -\o“a  is  the  Cauchy  pressure  and 

B  =  -J(dp/  dJ)  =  p(dp  /  d  p)  (7) 

is  the  bulk  modulus,  in  (6)  measured  at  J  =  1 . 
Properties  for  SiC  and  AlON  single  crystals  are 
listed  in  Table  1  [11-15], 

TABLE  1.  Elastic  properties  for  SiC  and  AlON  crystals. 


Property2 

Value  (SiC)b 

Value  (A10N)C 

Ci, 

5.01 

3.01 

C\2 

1.12 

1.55 

C44 

1.61 

1.74 

Q3 

0.52 

(=c12) 

C33 

5.49 

(=Cn) 

Cm 

-8.44 

-10.39 

Cm 

-8.73 

-7.13 

Cm 

-8.73 

-7.13 

C123 

-8.88 

-5.50 

C133 

-8.73 

-7.13 

C144 

0.15 

-0.81 

C155 

0.15 

-0.81 

C222 

-8.44 

-10.39 

C333 

-8.44 

-10.39 

C344 

0.15 

-0.81 

Bulk  modulus  B 

2.22 

2.04 

Shear  modulus  G 

1.94 

1.34 

Poisson’s  ratio  V 

0.16 

0.23 

Density  p0  [g/cm3] 

3.227 

3.714 

■“Elastic  constants  Cap  and  Capx  in  units  of  Mbar; 

Eq.  (6)  used  for  Capx. 

^References  [11,  13].  “References  [12,  14,  15]. 


A  cohesive  zone  model  addresses  intergranular 
fracture,  i.e.,  stress-induced  separation  between 
grains  in  a  given  polycrystal  [14,  16-20].  Each 
grain  in  a  polycrystal  is  treated  as  a  distinct  solid 
body  interacting  via  contact  with  surrounding 
grains  [21].  Prior  to  attainment  of  critical  traction 
tc,  contact  between  grains  is  rigid  (i.e.,  perfect 
bonding).  When  stresses  are  sufficient  such  that  the 
critical  traction  is  achieved  at  a  given  interface, 
separation  ensues,  with  traction  held  at  the  critical 
level  until  complete  separation  occurs  beyond 
distance  Sc.  Analogous  relationships  and  identical 
material  properties  are  used  for  mode  I  (normal) 
and  mode  II  (tangential)  fracture.  Bulk  elastic  and 
surface  properties  are  related  by  [17,  20] 

tcSc=K2(l-v2)/E  =  2r,  (8) 

where  K  is  the  fracture  toughness,  E  =  2G(l+v)  is 
Young’s  modulus,  and  F  is  surface  energy.  After 
complete  separation  occurs  in  either  normal  or 
tangential  fracture,  atomic  bonds  are  considered 
irreversibly  broken,  and  boundary  facets  interact 
subsequently  via  locally  frictionless  contact  with 
material  interpenetration  prohibited  [21].  Fracture 
properties  are  listed  in  Table  2  [22,  23]. 

TABLE  2.  Fracture  properties  for  SiC  and  AlON. 


Property 

Value  (SiC)a 

Value  (A10N)b 

to  [GPa] 

0.570 

0.306 

4  [pm] 

0.10 

0.06 

r  [J/m2] 

28.1 

9.0 

K  [MPa  m1/2] 

5.1 

2.5 

“References  [6,  22]  .  bReference  [23]. 


MICROSTRUCTURES 

In  the  absence  of  microstructures  from  actual 
test  specimens,  finite  element  meshes  of  synthetic 
microstructures  are  used.  Polyhedral  grain 
geometries  are  produced  using  a  Monte  Carlo  grain 
growth  algorithm  [24].  Volume  meshes  of 
tetrahedral  elements  are  created  from  stereo¬ 
lithographic  (STL)  files  of  surface  representations 
of  grains  comprising  a  given  polycrystal  [20]. 


TABLE  3.  Features  of  polycrystalline  microstructures. 


Feature 

Microstructure  I 

Microstructure  II 

Dimensions 

lxlxl  mm3 

lxlxl  mm3 

No.  grains 

50 

126 

Grain  size 

270  pm 

200  pm 

No.  elements 

1.59xl06 

1.13X106 
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Two  microstructures  are  considered,  with 
features  listed  in  Table  3.  Material  properties  of 
SiC  or  AlON  are  assigned  to  either  microstructure 
in  a  given  simulation.  In  this  way,  effects  of  elastic 
and  fracture  properties  on  the  overall  response  are 
studied  while  keeping  the  microstructure  (i.e., 
mesoscale  geometry)  fixed.  Loading  in  different 
directions  provides  an  indication  of  effects  of  grain 
structure  on  anisotropy.  Furthermore,  different  sets 
of  random  initial  lattice  orientations  are  assigned  to 
polycrystals  among  different  simulations,  enabling 
quantification  of  effects  of  anisotropic  elasticity. 

NUMERICAL  SIMULATIONS 

Lagrangian  finite  element  calculations  are 
performed  using  the  SIERRA  explicit  dynamics 
code  [21].  Considered  first  are  results  for  uniaxial 
strain  compression.  A  uniform  velocity  gradient  is 
assigned  as  an  initial  condition.  The  applied  strain 
rate  is  £  =  105/s.  Under  these  conditions,  the 
current  volume  V  of  the  specimen  is  related  to  its 
initial  volume  V0  via 

V  =  V0(1  -et).  (9) 

Figure  1  shows  axial  stress  in  a  SiC  specimen  at 
6%  compression. 


Figure  1.  Axial  Cauchy  stress  (positive  in  compression), 
microstructure  I,  SiC,  uniaxial  strain,  V /V0  =0.94  . 


In  uniaxial  strain  simulations,  one  face  of  the 
microstructure  is  assigned  a  constant  velocity  in  the 
compression  direction  while  lateral  faces  are 
restricted  from  expanding  or  contracting.  Evident 
are  stress  concentrations  near  grain  boundaries  and 
axial  splitting  cracks  typical  of  brittle  materials 
with  low  Poisson’s  ratio  [25].  Figure  2  shows 
average  axial  stress  E  computed  for  representative 
simulations  with  nonlinear  or  linear  elasticity  as 

Z  =  ±\\t^dA\,  (10) 

where  t'"]  is  the  component  of  traction  normal  to 
the  compressed  face  of  the  microstructure  with  area 
A.  Also  shown  for  comparison  are  experimental 
plate  impact  data  [26-29].  Linear  elasticity  gives 
smaller  normal  stress  than  corresponding  nonlinear 
models.  Figure  3  shows  average  shear  stress 

r  =  i|2-i-2-3|,  (11) 

with  { 27, ,  27,  J  the  {maximum,  minimum}  principle 
stresses  computed  similarly  to  (10).  Predictions  for 
SiC  are  generally  lower  than  experiment;  those  for 
AlON  are  generally  closer  to  experimental  data. 


v/v0 


v/v0 

Figure  2.  Average  axial  stress,  uniaxial  strain. 
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v/v0 


v/v0 

Figure  3.  Average  shear  stress,  uniaxial  strain. 


Unconfmed  compression  simulations  are  also 
performed  on  SiC  and  AlON  microstructures,  with 
initial  lattice  orientation  distributions  and  loading 
axes  varied  among  simulations.  These  and  all 
subsequent  results  incorporate  nonlinear  elasticity. 
Boundary  conditions  again  correspond  to  a  strain 
rate  of  £  =  105/s,  with  a  uniform  initial  velocity 
gradient  prescribed.  Peak  compressive  stresses  If, 
computed  as  in  (10),  are  compared  with  Kolsky  bar 
data  [30-32]  in  Table  4.  Again,  predictions  for 
AlON  are  closer  to  experimental  data;  however, 
loading  rates  in  experiments  (~103/s)  are  smaller 
than  those  applied  in  the  numerical  simulations. 


TABLE  4.  Unconfined  compressive  strength  2)  [GPa], 


SiC 

AlON 

Model 

4.2-4.9 

3. 1-3.8 

Experiment 

5.0-8. 2a 

3.0-4. 0b 

“References  [30,  31]. 

Reference 

[32], 

conditions  provide  a  shear  strain  rate  of  y  =  1 0"  /s  , 
with  a  uniform  initial  velocity  gradient  prescribed. 
In  some  simulations,  all  faces  are  free  to  expand 
laterally  to  accommodate  dilatation;  in  others,  the 
face  on  which  the  shear  traction  is  applied  is 
prohibited  from  expanding,  thereby  increasing  the 
confining  pressure  that  resists  dilatation  and 
impedes  shear  [2-4].  The  former  conditions  are 
labeled  “free”;  the  latter  are  labeled  “confined.” 
Figure  4  shows  the  shear  stress  component 
conjugate  to  the  applied  strain  rate  component  in 
representative  simulations.  Local  shear  stresses  are 
higher  for  confined  boundary  conditions;  similar 
trends  are  observed  for  average  shear  stresses. 


Figure  4.  Shear  stress,  microstructure  II,  AlON, 
y  =  0.03  ,  displacements  magnified  5x:  unconfined  shear 
loading  (upper  image),  confined  shear  loading  (lower). 


Shear  simulations  are  also  performed  on  SiC 
and  AlON  microstructures,  with  initial  lattice 
orientation  distributions  and  loading  planes  and 
directions  varied  among  simulations.  Boundary 


Failure  statistics  from  numerous  unconfmed 
shear  simulations  (24  for  SiC;  19  for  AlON)  are 
now  analyzed.  Average  shear  stress  T  is  computed 
analogously  to  (10)  for  a  given  polycrystal: 
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T  =  ^\\t{ndA\,  (12) 

where  l"'  is  the  tangential  component  of  traction 
applied  along  A.  Shear  strength  is  the  (peak)  value 
of  rat  which  dz/dy=  0.  Table  5  shows  normal 
(i.e.,  standard  Gaussian)  statistics  computed  from 
discrete  values  of  r  predicted  for  each  material 
from  simulations.  The  mean  value  of  r  is  denoted 
by  f ,  standard  deviation  (square  root  of  the 
variance)  by  S.  Mean  strength  is  lower  in  AlON 
than  SiC,  presumably  a  result  of  its  lower  cohesive 
strength,  surface  energy,  and  elastic  stiffness. 
Standard  deviation  is  larger  in  AlON  than  SiC, 
possibly  a  result  of  greater  anisotropy  and  smaller 
fracture  strength  and  energy  that  render  AlON 
more  sensitive  to  cracking  at  grain  boundaries  or 
triple  junctions  favorably  oriented  for  fracture. 


TABLE  5.  Normal  statistics:  unconfined  shear  strength. 


SiC 

AlON 

f  [GPa] 

5  [GPa] 

T  [GPa] 

5  [GPa] 

Model 

0.511 

0.013 

0.333 

0.027 

Table  6  lists  Weibull  statistics  for  values  of  r 
predicted  from  simulations;  shown  for  comparison 
are  statistics  from  static  flexure  experiments  on 
much  larger  specimens  [5,  34].  A  cumulative 
Weibull  distribution  predicts  probability  / of  failure 
at  or  below  a  given  shear  stress  T  as 

/(r)  =  l-exp[-(r/ron,  (13) 

with  nominal  strength  r(J  and  Weibull  modulus  m. 
Larger  m  corresponds  to  lower  variability.  Weibull 
parameters  are  computed  from  discrete  data  points 
by  linear  fitting  to  In  [ln[l/(l-/)]}  versus  lnr, 
e.g.,  following  [33].  Predictions  for  SiC  exhibit 
lower  variability  and  greater  nominal  strength  than 
those  for  AlON,  in  qualitative  agreement  with 
experiments.  However,  because  the  mesoscale 
models  do  not  incorporate  initial  flaws  or  spatially 
variable  fracture  properties,  the  present  results  may 
under-predict  variability. 


TABLE  6.  Weibull  statistics:  unconfined  shear  strength. 


SiC 

r0  [GPa] 

m 

AlON 
r0  [GPa] 

m 

Model 

0.52 

45.2 

0.35 

14.9 

Experiment 

0.37- 

0.62a 

4.9- 

26.6a 

0.29- 

0.81b 

2.9- 

26.3b 

'‘Reference  [5]  bReference  [34] 


CONCLUSIONS 

Mesoscale  models  incorporating  nonlinear 
anisotropic  elasticity  and  cohesive  fracture  have 
been  used  to  study  dynamic  compression  and  shear 
of  SiC  and  AlON  ceramics.  Average  stresses 
predicted  for  compression  are  in  reasonable 
agreement  with  experiments  on  much  larger 
samples  of  material;  agreement  between  model  and 
experiment  is  generally  closer  for  AlON  than  SiC. 
Predictions  for  dynamic  unconfined  shear  strength 
demonstrate  greater  variability  for  AlON  than  SiC, 
in  qualitative  agreement  with  experimental  static 
fracture  statistics. 
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